Composite to tilted vortex lattice transition in Bi 2 Sr 2 CaCu 2 08+5 in oblique fields 
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Precision measurements of the vortex phase diagram in single crystals of the layered superconduc- 
tor Bi2Sr2CaCu2 0g+i in oblique magnetic fields confirm the existence of a second phase transition, 
in addition to the usual first order vortex lattice melting line H m (T). The transition has a strong 
first order character, is accompanied by strong hysteresis, and intersects the melting line in a tri- 
critical point (H^, Jf| r ). Its field dependence and the changing character of the melting line at the 
tricritical point strongly suggest that the ground state for magnetic fields closely aligned with the 
superconducting layers is a lattice of uniformly tilted vortex lines. 

PACS numbers: 74.25.Qt, 74.25.Op, 74.25.Dw 



The first order "vortex melting" transition from a solid 
(phase-ordered) state to a liquid state with only short 
range correlations is the main feature of the phase dia- 
gram of vortex lines in clean, layered high-temperature 
superconductors [l[ . The application of a small field com- 
ponent 7?H, parallel to the superconducting layers, leads 
to a lattice of tilted vortex lines that melts in a similar 
fashion However, in the more anisotropic (layered) 
compounds such as Bi2Sr2CaCu20s+5, the depression of 
the perpendicular component of the melting field by 
larger parallel fields was interpreted as the consequence 
of the decomposition of the tilted vortex lattice into a 
combined lattice structure of Josephson Vortices (JVs) 
and Abrikosov-type pancake vortices (PVs) Q ■ For very 
small field components H perpendicular to the layers, 
chain structures Q arising from the attractive interac- 
tion of PVs with JVs were directly visualized by Bitter 
decoration scanning Hall-probe @,@] and magneto- 

optical techniques @, S] ■ At higher H 1 - ~ , the con- 
tribution of the JVs to the free energy of the pancake 
vortex crystal results in the almost linear dep ression of 
as function of the parallel field [1 EE EI El ■ This 
behavior in moderate ff" stops at a temperature depen- 
dent characteristic field H c \-. Even though melting is still 
observed above H cr , the variation of with increasing 
if II becomes much weaker 
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Several controversial 
interpretations of this ch ang ing behavior were proposed, 
such as layer decoupling [ll |j_a commensurate transition 
[3], and a matching effect 3). 

In this Letter we focus on the high-temperature por- 
tion of the vortex phase diagram in single crystalline 
Bi 2 Sr 2 CaCu208+5 in oblique fields, which can be estab- 
lished precisely using the well-defined discontinuity of the 
vortex density at the melting transition. We show that 
(H^, Her) corresponds to a tricritical point in the vortex 
lattice phase diagram, where the melting crosses a novel 
transition from a composite lattice at low parallel fields, 



to another tilted lattice structure at high H" . The exper- 
imental observation of large hysteresis suggests that this 
transition is strongly first order, consistent with recent 
predictions [ToT ] . The identification of the vortex ground 
state at high parallel field as a tilted lattice structure re- 
solves the open problem of the apparent anisotropy fac- 
tor Jeff, and allows one to determine the enhancement 
of by magnetic coupling. We find the temperature 
dependence of y e / / to be consistent with previous obser- 
vations 3, 17 1 and in quantitative agreement with the 
proposed model. 

Experiments were performed on rectangular samples 
cut from Bi2Sr2CaCu208+<5 single crystals with different 
oxygen content The c-axis component of the local 

magnetic induction B (r) was measured by micro-Hall 
sensors placed on the central part of the sample. The 
2D electron gas Hall sensors were fabricated in GaAlAs 
heterostructures and had 8x8 fim 2 active area. Re- 
sults are presented in Fig. [Ha) as the local magnetiza- 
tion H^- = B^ — . The local dc magnetization of all 
crystals shows a sharp discontinuity, AB 1 - , at the vortex 
melting transition, that was tracked as function of if" 
at various fixed temperatures. The angle 9 between the 
magnetic field and the crystalline c-axis was computer- 
controlled with 0.001° resolution, while the field magni- 
tude could be swept up to 1 T using an electromagnet. 
Two types of magnetization loops were measured. In the 
first, the magnetization is traced as function of the c-axis 
field at constant if"; in the second, the magnetization is 
measured as function of if" at constant H ± . 

While the discontinuity in the dc magnetization gives 
a clear identification of the melting field, another method 
[lit , in which the magnitude B(f, T) of the periodic part 
of the induction above the sample is measured at the fre- 
quency / of an ac ripple field applied perpendicularly to 
the sample plane, is more convenient and precise. The 
ac response is represented as the transmittivity T', i.e. 
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FIG. 1: (a) dc local magnetization loops recorded on an as- 
grown Bi2Sr2CaCu208+i single crystal (T c = 88 K) at 75 K, 
as function of the magnetic field component H perpendicu- 
lar to the superconducting layers, with the in-plane field Hu 
held constant. The inset shows a magnified view of the dis- 
continuity at the vortex melting transition in H" = 0. (b) the 
in-phase (screening) component of the ac response of the same 
crystal, recorded under the same conditions, with an ac mag- 
netic field of amplitude h ac = 0.8 Oe and frequency / = 11 Hz 
applied along the c-axis. The melting transition shows up as a 
paramagnetic peak (see inset), the transition from combined 
to tilted vortex lattice is indicated by arrows. The ac response 
is plotted as the transmittivity T 1 = [B'(f,T) - B(f,T < 
T c )]/[B(f,T > T e ) - B(f,T< T )). 



the in-phase component B'(f, T), normalized by the am- 
plitude h ac of the ac ripple [2l| . The steplike feature in 
the dc magnetization loop at translates to a param- 
agnetic peak in the ac response, shown in Fig. [TJb) . 
The magnitude of this peak depends on the ratio of AB 
to h ac . The peak position is independent of both the am- 
plitude and frequency of the ac ripple. In the explored 
temperature range (above 50 K) and at low frequency 
(below 27 Hz), a true paramagnetic signal is measured. 
At higher frequencies or lower temperatures, flux pinning 
results in the partial shielding of the ac field [20(. Nev- 
ertheless, a peak-like feature persists at melting. 

Figure HJa) shows that at T > 50 K, the application 
of even a small magnetic field component parallel to the 
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FIG. 2: Transmittivity T 1 of the same crystal as in Fig. [T] 
recorded at T — 70 K and constant ff x = 58 Oe, as a function 
of if" for various frequencies / and amplitudes h ac of the 
ac ripple field. A marked hysteresis of the ac screening is 
observed. This hysteresis disappears when h ac is increased. 



layers results in the drastic suppression of magnetic irre- 
versibility. This is expected when the geometric barrier 
is at the origin of flux pinning 2^, 23|. Simultaneously, 
H^ is depressed linearly with increasing £fH. However, 
at a well-defined value Her, the dependence of H^ on 



in-plane field changes to a much slower, quadratic be- 
havior that very well fits the anisotropic London model, 



H m {6) = ff m0 /(cos 2 6» + sm 
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The characteristic field H m Q and the effective anisotropy 
parameter j e f f will be defined below. 

In Fig.[ljb), another feature in the in-phase component 
of the ac response can be distinguished, at perpendicular 
fields H somewhat smaller than H^. This feature is 
brought out much more clearly in sweeps of the parallel 
field, shown in Fig. [2j There is an abrupt jump from 
lower to higher values of T' on increasing H " , that only 
appears for parallel fields J?" < i/J r . The position of the 
jump does not depend on ac frequency. At low amplitude 
of the ac field, a pronounced hysteresis of T' is observed; 
this disappears if h ac is sufficiently increased. 

The transmittivity T" is simply related to the mag- 
nitude of the shielding current flowing in the sample in 
response to the applied ac magnetic field, a higher T" cor- 
responding to a smaller current and less screening 21 1. 
In the present case, dc magnetization loops point to the 
geometrical barrier 22j as the main source of screening. 
However, increasing the ac field frequency reduces the 
role of thermally activated depinning of vortices in the 
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FIG. 3: Two field-component vortex-lattice phase diagram, 
with the first-order melting transition (o) determined 
from the paramagnetic peak in T , and the first-order tran- 
sition to the tilted PV lattice at Bt (•), determined from 
the "glitch" in T' (Fig. [Jb). The dashed line is a fit to the 
composite-to-tilted lattice transition, Eq. Q with C = 0.030; 
the continuous line is a fit of the high-field portion of the vor- 
tex lattice melting line to Eq. {TJ. 



crystal bulk; as a consequence, a bulk screening current 
due to vortex pinning emerges [20| . At the frequencies 
of Fig. [2j the step in T" is due to a discontinuous change 
of the magnitude of this bulk current at the well-defined 
in-plane field, ifH = H^ t . The location of does not 
depend on the frequency and h ac which indicates a vortex 
phase transition in the bulk, from a low ff"-phase with 
higher pinning, to a high ij"-phase with lower pinning. 
The hysteresis of the screening current indicates it to be 
first order. 

In Fig. [3] we collect, for T = 75 K, the positions of 
the two first order transitions in a plot of H versus 
H". The usual melting field of the vortices, de- 
duced from the paramagnetic peak in the transmittiv- 
ity, shows the well-known linear decrease as function of 
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13], up to the field component Her- The 
field H^. of the first order transition revealed by the (ir- 
reversible) transmittivity rapidly increases with iJy and 

crosses the melting line at Her in a tricritical point. The 
same scenario is observed at all explored temperatures 
(T > 50 K), with temperature dependent values of ifj r . 
The anisotropy factor 7 e //, extracted from the London 
model fit to the high-iP part of the melting line, de- 
pends on temperature as well as on the oxygen content 
of the Bi 2 Sr 2 CaCu208+5 crystals; it is depicted in FiglU 
The \ow-H " portion of the phase diagram with almost 
linear H^(H^) dependence has been interpreted as the 
region of crossing vortex lattices of JVs and PV stacks Q • 
A more accurate analysis shows that other "composite- 
lattice" configurations compete for the ground state in 
the parameter range of I^S^CaC^Os+a- These are 
the soliton lattice [25J , as well as the set of combined lat- 



tices composed of regularly spaced rows of tilted pancake 
stacks, separated by M rows of pancake stacks aligned 
with the c-axis. The latter type of lattice becomes fa- 
vorable at smaller anisotropies and larger H". Moreover, 
if H" is sufficiently large and the material anisotropy is 
not extremely high, a simple tilted lattice (M = 0) turns 
out to be the most favorable configuration. We inter- 
pret the experimentally observed transition as that from 
a composite to such a uniformly tilted lattice. A simple 
estimate for the in-plane field at which this transition is 
expected, B\ t , can be obtained by comparing the ground 
state energies of the simplest (M = 1) composite lattice 
and of the uniformly tilted lattice, giving [15| 



B-^o/ln 



1.55-/B -L$ 



-,1/2 



sB 



(2) 



Here <2>o is the flux quantum, A Q & is the a6-plane pen- 
etration depth, 7 is the penetration depth ratio A c /A a &, 
and s is the layer spacing. Eq. © gives a very good fit to 
the experimental transition line, as illustrated in Fig. [3] 

The anisotropic three-dimensional behavior (pj of 
for large in-plane field > Bar strongly supports this 
interpretation. The H^(H" )-dependence is the direct 
consequence of the vanishing contribution of the mag- 
netic interaction between PVs to the vortex tilt stiffness 
in a highly inclined tilted vortex structure. The angular 
dependence of the melting field can be derived using a 
scaling transformation of coordinates, z = 7 2 / 3 z; r±_ = 
-y _1 / 3 rj_, which reduces the larger part of the free energy 
to an isotropic form j^j • In scaled coordinates the mag- 
netic field is given by B = Bj 2 / 3 (cos 2 6 + sin 2 fl/7 2 ) 1 ^ 2 , 
while the tilt angle tan (9 = tan O/'y. The Josephson tilt 
energy of a deformed vortex line (PV stack) in scaled 
coordinates, 
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is determined by the effective line tension i\(ki) = 
Sq In (l/r cn ^ki V valid when the wave vector along the 

line direction, hi, is much larger than the vortex lattice 
zone boundary vector. Here, 6u(ki) is the Fourier trans- 
form of the line deformation, and io = £o7 -2 ^ 3 with Eq = 
<I>g/(47rA a 6) 2 . For near-perpendicular fields (ffel) the 
core cut-off distance f cu t is determined by the so-called 
thermal vortex wandering length, f cut w (u 2 n+1 ) 1 / 2 = 
((Gn+i — Un,) 2 ) 1|<2 , where u„ is the position of the PV 
vortex in layer n [26j | . For a tilted vortex line, u„ in+ i = 
s tan# + (5u n ,„ + i consists of the average displacement as 
well as random (thermal) fluctuations meaning that the 
core cut-off f 2 ut ps s 2 tan 2 6 + ((<5u„.„+i) 2 )- The melt- 
ing temperature is given by T m = A^/ii(l / a)ioa, with 
a = (Qa/B) 1 / 2 and ArsO.1 [27[. Returning to real coor- 
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FIG. 4: Temperature dependence of the apparent anisotropy 
y e ff extracted from the fits of the melting line to Eq. fl}. 
The drawn line shows a fit to Eq. (j4j) with intrinsic 7 = 500. 



a strong first order character and intersects the usual 
first order vortex lattice melting line at a tricritical point 
[H^{T),H l } t (T)}. For fields parallel to the superconduct- 
ing layers if" < H cr the melting line shows the signature 
of a composite lattice. For if 1 1 > Her, the melting line 
is fully consistent with that of a uniformly tilted lattice 
of PV stacks. We thus propose that the new first or- 
der transition takes place between the combined and the 
tilted vortex lattice. For low in-plane fields, the com- 
bined vortex lattice is stabilized by the magnetic interac- 
tion between PV's in the same stack (vortex line). The 
enhancement of the melting line in the combined lattice 
regime is due to the contribution of this magnetic inter- 
action to the vortex line tilt stiffness. 



dinates, we obtain 

T 2 _ 42/ s 2] ( CjB sc {6)/B \ B sc (6) 

where the numerical constant Cj ps 5 can be esti- 
mated within the self-consistent harmonic approxima- 
tion, B sc (9) = ($ /7 2 s 2 )/ Vcos 2 9 + 7-2 sin 2 0, and r 2 , = 
<(5un,„+i) 2 )/(7s) 2 « 2A[$ /s 2 7 2 S m (0)] 1 /2. Note that 
the angular-dependent core cutoff introduces an addi- 
tional angular dependence of melting field: T m no longer 
depends only on the ratio B sc (9) / B. In particular, a new 
angular scale appears given by tan# = 77*0. In the exper- 
imental angular range tanf? <C 7,77-0, we recover Eq. ([T]) 
with the apparent anisotropy 

^ / KViUoh 2 * 2 /^ \ 1/2 

Jeff ~ 7 1 + — F =T • (4) 

\ In [68 v /$o/( J B„ i (0)7 2 s 2 )J J 

We note several key points. First, the effective anisotropy 
Jeff is manifestly smaller than the intrinsic 7. It in- 
creases with temperature, and is in excellent agreement 
with the experimental data of FigJH strongly suggesting 
that the modified core cut-off length originating from the 
tilting of the PV stacks determines the behavior of the 
melting line at high parallel fields. Very similar behavior 
has been observed in YBa2Cu307_^ [16|. Next, the pref- 
actor H^o = B m (0)/fi O m Eq. fl]) is to be interpreted 
as the hypothetical vortex melting field H m .j(6 = 0) in 
the absence of the magnetic coupling between PVs. The 
difference AH mag = H^{6 = 0) - H m ,j(0) « 0.15^ 
between the real (experimental) melting field and this 
prefactor represents the (remarkably modest) enhance- 
ment of the melting field due to magnetic coupling. 

Summarizing, we have established the existence of 
phase transition of the vortex lattice in single crystalline 
Bi2Sr2CaCu208+<5 in oblique fields. The transition has 
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